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Abstract 

We study descendants of inhomogeneous vertex models with boundary reflections 
when the spin-spin scattering is assumed to be quasi-classical. This corresponds 
to consider certain power expansion of the boundary- Yang-Baxter equation (or re- 
flection equation). As final product, integrable su(2)-spin chains interacting with 
a long range with XXZ anisotropy are obtained. The spin-spin couplings are non 
uniform, and a non uniform tunable external magnetic field is applied; the latter 
can be obtained when the boundary conditions are assumed to be quasi-classical 
as well. The exact spectrum is achieved by algebraic Bethe ansatz. Having real- 
ized the su(2) operators in terms of fermions, the class of models we found turns 
out to describe confined fermions with pairing force interactions. The class of mod- 
els presented in this paper is a one-parameter extension of certain Hamiltonians 
constructed previously. Extensions to su(n)-spin open chains are discussed. 



1 Introduction 



Integrable vertex models (VM) in two dimensional classical statistical mechan- 
ics are the common seed of many relevant exactly-solved quantum models in 
one dimension [1,2]. Famous examples are the XXX, XXZ, and XYZ Heisen- 
berg chains that find more and more applications in contemporary physics. 
The key toward this powerful synthesis is to notice that the "scattering" of 
the degrees of freedom of both the VM and the spin chains is described by 
the same matrix. The Quantum Inverse Scattering Method (QISM) exploits 
this fact systematically [3] . The method relies on the observation that transfer 
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matrices t = Tr(T) span a one-parameter family of commuting operators if a 
(scattering) matrix R exists such that T, R satisfy the celebrated Quantum 
Yang-Baxter equation. The equivalence between VM and Heisenberg chains 
consists in the fact that these models have the same i?-matrix. Due to the 
property of the scattering R(u, v) = R{u — v) , Viz, v e C the integrabil- 
ity of VM is preserved if disorder is added at each lattice site such that the 
scattering "wave momenta" iz, v result to be shifted arbitrarily. In this case, 
however, it is difficult to extract a Hamiltonian. A route to simplify the prob- 
lem is to resort the so called "quasi-classical" limit of the QISM. The term 
"quasi-classical" here indicates that the scattering between the degrees of free- 
dom of the model is assumed to be quasi-classical. Quantitatively, this means 
that a parameter 77 does exist such that the scattering matrix is of the form 
R(u) oc 11 <S> 11 + i]r(u) in the limit i] — > (rj plays the role of K). The quantity 
r(iz) fulfills the classical Yang-Baxter equation (that is a restatement of the 
Jacobi identity for the Poisson brackets of suitable action-angle variables). It 
is worthwhile to mention, however, that the systems obtained by this quasi- 
classical expansion consist of quantum spins (by no means quasi-classical). 
The quasi-classical expansion of the transfer matrix of disordered VM (in the 
lowest spin representation) is non trivial and it produces the Gaudin's magnet 
Hamiltonians[4,5] containing a long range spin interaction (in contrast with 
the range of the Heisenberg chains which involves nearest neighbour spins). 
A richer variety of integrable models by QISM comes from imposing non trivial 
boundary conditions different from the periodic ones. Twisted boundary con- 
ditions, for example, imposed to the six vertex model [6] produce the Gaudin 
magnet in a non-uniform local magnetic field, which is very important for 
physical applications. In fact having realized the (pseudo)spin algebra in terms 
of fermions the XXX Gaudin Hamiltonians in a non uniform magnetic field 
are the constants of the motion of the BCS model [7] that describes pairs of 
electrons (in time reversed states) interacting with a long range uniform pair- 
ing coupling. The exact solution of the BCS model was found long ago by 
Richardson [8] and rediscovered recently. In particular it was used to study 
small metallic grains [9,10]; the picture was merged in the scenario of QISM 
in the Ref. [11]. Connections with WZNW models in field theory have been 
deeply investigated [12] based on the relation between solution of KZ equation 
and Gaudin model found in Refs. [13,14]. The class of pairing Hamiltonians 
was generalized by investigating the quasi-classical expansion of the disor- 
dered twisted six vertex model with XXZ i?-matrix. In terms of fermions 
this class of Hamiltonians represents interacting electron pairs with certain 
non-uniform long-range coupling strengths [15,16,17,18]. Twisted rings can be 
cut to open chains and loops include two reflections at the boundaries. The 
possibility to include such reflections in integrable theory was founded and 
systematically investigated by Sklyanin [19]. The quasi-classical limit of the 
disordered six vertex model with boundaries was investigated first by one of 
the authors [20,21]. This led to a model where the spin couplings contain an 
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additional parameter with respect to the original Gaudin magnet, and in a 
vanishing external magnetic field (see Eq. (27) and the relative discussion be- 
low of it). In the present work, we proceed along this line. We still consider an 
inhomogeneous six-vertex model with boundary reflection, following closely 
Ref. [20]. By properly choosing the reflection parameters, we introduce an 
external non-uniform magnetic field of tunable strength in the Hamiltonian. 
The trick consists in the assumption that also the boundary conditions have a 
quasi-classical expansion (see Eq. (25) and Sec. 3.2). At best of our knowledge, 
this idea is pursued for the first time in the present paper. In the following 
we summarize the main results obtained in the bulk of the paper. The class 
of spin-Sj (j = 1, . . . , N) models that we find has Hamiltonian of the form 

h = £2^5; - £ [jgsisi + (J jk s+sz + h.c.)] . (i) 

We agree that the latin indices j, k will run from 1 to N, where N is the 
number of spins. The operators S ± , S z are su(2) operators. The couplings 
are 

Jjk = (cosh(2pzj) + cosh.(2pz k ) — 2 cos (2pt)^j 
Jjk = Ijk (sinh[p(z,- - it)} sinh[p(z k + it)]) 

j = j(s z ) hj - hk 

jk cosh(2pZj) — cosh(2pz k ) 

J(S z ) = j(l-JS z )~ 1 (2) 

where S z is the total ^-component of the spin. The quantities hj,Zj are two 
arbitrary sets of real parameters; t is also a real arbitrary parameter and it 
directly comes from the boundary terms (see Eq. (16) with £ = it); finally, p 
can be 1, i or can be tending to zero corresponding to hyperbolic, trigonomet- 
ric, and rational couplings, respectively. 

The eigenstates in the sector with total ^-component of the spin S z = J2j Sj — 

M, are 

M 

m=\{S-{e a )\H) , (3) 



a=l 



where 



and 



S~ (u) = V cosh b( M + z i + 2it )\ ~ cosh[p(M - Zj )] ^ 

[) y cosh(2pu) - cosh(2pZj) j {) 



N 

\H) = ® 
The corresponding eigenvalues are 



S z = S, 



£ = E 2 Vi (5) 
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where 



T] = s 3 (i - j(s') £ s k 1 —^ + j(s') £ Lj^sl] ? (6) 



/■ r , X j a •' J A,, 



where we defined 

cosh(2pZj) — cos(2pt) , — ^— - — = cosh(2pe a ) — cos(2pt) (7) 



— - — - = cosh(2pZj) — cos(2pt) , ~*~ ^ a 



The rapidities \ a satisfy, in the sector having total ^-component of the spin 
S z , the Bethe equations: 

1 <?■ 

E — E + 

i / l + JQS'Xi + g') l- J(g')(i + g») \ = u 

2 J(S*) \ 1 + X a 1-X a J 

We point out that the dependence on the reflection parameters comes only in 
the couplings and eigenvectors (Eqs. (2) (3)), while the eigenvalues depend on 
t only implicitely (through Eq.(7)). 
The rational limit of the models is recovered for p — > 0. 

For t — and pt = tt/2, the models reduce to the ones that we presented in 
Ref. [15] (see section 3.3). Thus the class of models we discuss in the present 
paper is a one-parameter extension of the former class. 

Using the fermionic realizations of su(2) the Hamiltonian (1) can be rephrased 
to describe confined fermions interacting with pairing and exchange forces (see 
Eq. (48)). 

The paper is organized as follows. In the next section we summarize the main 
ingredients of the inverse scattering of VM with boundaries. In section III 
we construct the integrable models we deal with together with their exact 
solution. In section IV we use the fermionic realization of the su(2) algebra to 
rewrite the Hamiltonians in a second quantized form. Section V is devoted to 
final remarks. In appendix A we review basic properties of VM. In appendix B 
we prove the integrability of a class of models when a more general (off- 
diagonal) reflection at the boundary is applied (see Eqs. (B.l) - (B.4)). We 
also discuss a generalization to su(n) case in appendix C. 



2 Integrable boundary conditions 



In this section, we review how the QISM is applied to VM, in order to obtain 
a family of commuting transfer matrices. VM describe a system of interacting 
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classical objects on a two dimensional lattice. As described in appendix A, 
the partition function of the system can be written as Z = Tr{£(l) ■ • -t(K)}, 
where t(i) are operators in some appropriate many-body linear space (in the 
sense that it is the direct product of N elementary linear spaces). The VM is 
exactly solvable if = 0. Usually, it is assumed that the dependence 

on the i-th row of the lattice comes through a parameter Ui, which takes 
values on some domain of the complex plane. Then the requirement for exact 
solvability becomes t(u),t(v) = 0, Vw, v belonging to the domain. 

The QISM provides a way of constructing classes of commuting operators 
t(u), finding their eigenvalues and their common eigenstates, and extracting 
Hamiltonians whose integrals of motions are t(u). The QISM is a procedure 
which starts from the i?-matrix and from the Lax operator to yield the transfer 
matrix t(u). From the transfer matrix, a class of Hamiltonians can be extracted 
in various ways, to be depicted below. The QISM has a built-in Algebraic 
Bethe Ansatz (ABA) which provides the diagonalization of the t(u), and hence 
of the Hamiltonian. 



The XX Z iZ-matrix is 



where 



R(u, v) 



U(u,v) \ 
b(u, v) c(u, v) 
c(u, v) b(u, v) 

\ a(u,v)J 



(9) 



a(u,v) = sinh \p(u— v+r))]/p , b(u,v) = sinh \p(u—v)]/p , c(u,v) = sinh.(pr])/p. 
It is connected to the i?-matrix defined for VM by R(u, v) = VuR^l, 2), where 

(\ 0^ 
10 
10 
^0 l) 

is the permutation operator, and it is assumed that the dependence of R upon 
the rows comes through a parameter assigned to each row. 
The corresponding Lax operators are 



V 12 = 



I ( smh[p(u + r]S*)} sinh(p?7) 
LAu) — - I 



P 



smh(prj) Sj~ smh[p(u — r/Sj)} 



(10) 
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Here p is the anisotropy parameter, in the sense that, when p ^ — in which 
case one can put either p = 1 or p = i — the QISM yields a Hamiltonian with 
XXZ-type couplings, while in the limit p — > 0, the hyperbolic/trigonometric 
functions reduce to rational ones, and the QISM generates a Hamiltonian 
having XXX couplings; rj, instead, is the so-called quantum parameter which 
plays the role of h; as we shall later see, it gives the degree of deformation of 
the classical algebra su(2) into the quantum algebra su q (2). We remark that 
the terminology is somehow misleading: since we associate the algebra su(2) 
with spins, realized either by true spins or by pairs of time-reversed electrons, 
in the limit r\ — > we obtain genuine quantum Hamiltonians. 

The Lax operators act on the auxiliary two-dimensional vector space V, and 
on the quantum space H n . They obey the fundamental Yang-Baxter relation 
Eq.(A.l), which in terms of the R- matrix now reads 

12 2 1 

R(u - v) Lj(u - Zj) Lj (v - Zj) =Lj (v - z 5 ) Lj (u - Zj) R(u - v) , (11) 

Due to the additive property of the i?-matrix R(u, v) = R(u — v), parameters 
Zj taking into account on-site disorder through the lattice can be introduced. 

1 2 

As customary Lj(u) = Lj(u)®l, and Lj(u) = 11 (g)Lj(u); the external product 
is meant between two copies of the space V, while the multiplication of the 
elements of Lj, which are operators on Tij, is an internal product. The relation 
(11) is actually obeyed only for 1/2 spins, i. e. for dim(Ti,j) = 2. The order 
of the representation (that is the dimension of Hj) can be extended to larger 
values, keeping the dimension of V fixed to 2; however, one has to renounce 
to the algebra su(2), and introduce rather the quantum algebra su q (2), which 
ensures that the relation (11) is obeyed whatever is the representation of the 
algebra. The parameter q is related to the parameters p and rj by q = exp{prj). 
The commutation rules are 

oz oil _ c± . \c+ c-1 _ smli^^D 
^J-^ . [V;J- sinh^) ' [LZ) 

In the quasi-classical limit i] — > 0, or in the isotropic limit p — > 0, su q {2) 
reduces to su{2). 

Next, we consider the monodromy matrix T(u) = Li{u — Z\) ■ • ■ L N (u — z^). 
We have an internal product over V and an external one over Hj and Hj'] 
thus T(u) is an operator over V ®H\® ■ ■ ■ ® Hn- It has the form 



T(u) 



A{u) B{u) 
C{u) D(u) i 

with A, B, C, D operators over H. = ®j Hj. 



The local relation (11), and the ultra-locality property, [//•*(«), L c k d (v ) 



= 
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for j 7^ k, imply that T(u) fulfills the global Yang-Baxter equation 

R(u - v) T(u) T (v) =T (v) T (it) R(u - v) , (13) 

In the case of periodic boundary conditions (on the auxiliary matrix space), the 
quantities Tr{ T(u)} = A(u) + D(u), where the trace is on the auxiliary space 
V, generate a one parameter commuting family of operators which underlies 
an integrable model. 

Remarkably, the property of integrability is preserved for a wider class of non- 
trivial boundary conditions. Integrable boundary conditions are introduced by 
the so-called "boundary ^-matrices" satisfying the reflection equations [22,19] 

R(u - v) K- (u) R(u + v) K- (v) = 

K- (v)R(u + v) K- (u)R(u - v) (14) 



R{-u + v)K + {u)R{-u -v- 2r])K + (v) 



K + (v)R(-u -v- 2t])K + {u)R{-u + v) 



(15) 



Among the solutions of the reflection equations, we consider the diagonal 
ones, which yield Hamiltonians preserving the total ^-component of the spin 
S z . They depend upon the free parameters £± 

K+(u) = K(u + ri,S + ) 



where 



1 / sinh \p(u + £)] 



P 







sinh [p(u - 0] j 



The family of commuting transfer matrices is 

t(u) = Tr{K + {u)U{u)} 
( 



A(u) B{u) 
\C{u) V{u) 



where U(u) = T(u)K^u)T-\-u) 
monodromy matrix is defined as [3] 

T-\-u) = <r y T\-u + n)<r v det^ 1 T(-M + 77/2) , 



(16) 



(17) 



The inverse of the 



where a y is the Pauli matrix in the representation where cr z is diagonal, and 
det q T(u) = A(u - rj/2)D(u + rj/2) - C(u - r]/2)B(u + rj/2) is the quantum 
determinant, which is a sw(2)-number. 
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The eigenvectors of i(u), in the sector with S z = Y.j Sj ~ are given by [19] 



M 



n *m \h) , 

a=l 



where 



\H) = 



S* = S, 



is the pseudo-vacuum state having all maximum Sj eigenvalues. The eigenval- 
ues are 

= sinh[2p(n + 7y)] / _ 2 x x 

sinh[p(2w + 77)] \ L v ^ 
sinh[p(w - e Q - 77)] sinh[p(u + e a )] 



a(u)d(—u + rj) + 



\ sinh[p(u + e a + 77)] sinh[p(w - e a )\_ 

sinh(2p-u) / \ 
- sinhb(2 M + r ? )] i COSh [2p(^ - * + v)} ~ cosh(2p5)) 



(18) 



x 



X 



■ M 

n 

.0=1 



sinh[p(w — e Q + 77)] sinh[p(-u + e a + 2rj)] 
sinh[j9(-u + e a + 77)] sinh[p(u — e a )] 



a(— if + r))d(u) , 



where we put £ + + = 2<r, £ + — £_ = 25 and a(-u) = llj=i sinh[p(-u — 
Zj + r)Sj)]/p, d(u) = EljLi sinh[p(u — ^ — r)Sj)]/p. The e a satisfy the Bethe 
equations: 



cosh[2p(e a + cr)] — cosh(2p5) 



n 



sinh[p(e a - ep - 77)] sinh[p(e a + ep) 



cosh[2p(e a - a + 77)] - cosh(2p5) ^ sinh[p(e a - ep + 77)] sinh[p(e a + ep + 2rj)\ 
sinh[p(e Q - Zj - rjSj)] sinh[p(e Q + Zj - r)(Sj + 1))] 



n 



- 1 sinh[p(e a - Zj + rjSj)] sinh[p(e a + Zj + r](Sj - 1))] 



(19) 



The final step to obtain integrable models from the procedure above is to 
observe that transfer matrices can be used as generating functional of Hamil- 
tonians. A possibility is 

(20) 

u=u c 

In the homogeneous = Vj, Heisenberg Hamiltonians with nearest- 

neighbour interaction are obtained for u c — 0. The presence of disorder makes 
the application of Eq.(20) quite difficult. A particular value of u c for which 
the calculations can be done is u c — > 00; in this case the interaction in the 
Hamiltonian is long range [23]. Another way to face the problem is to resort to 
the quasi-classical expansion. The trick consists in obtaining a set of commut- 
ing operators as coefficients of the power-77 expansion of t(u) — from which 
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a Hamiltonian turns out can be built as a polynomial. The quasi-classical 
expansion of the i?-matrix and Lax operators gives 

l + rjl < J 1) (u)+rj 2 lf\u) + 0(rj 3 )] , 
R(u) = 1 <g> l + 77 r(u) + 0{t] 2 ) , 

where 11 is the (2 x 2) identity and r(u) reads 

N 

1 

1 

cosh(pu) y 

Unfortunately, if one wants to extend the results to spins higher than 1/2, one 
has to increase the dimension of the auxiliary space accordingly [24] . There is, 
though, the remarkable exception of isotropic models, i.e. the ones obtained in 
the p — >■ limit, for which the quantum algebra reduces to su(2). The XXX 
model with higher spin was introduced in Refs. [24,25]. A central point of the 
approach is that, the quasi-classical limit 77 — ^ reduces the quantum algebra 
to ordinary su(2), whatever is the dimension of the quantum space Tij. This 
implies that the operators thus found are realized through spin operators 
Sf, and that they commute with each other for arbitrary spins. 



3 Quasi-classical expansion of vertex models with boundary reflec- 
tions 

In this section we investigate systematically the power ^-expansion of the 
transfer matrix of disordered vertex models with boundary. As final product 
of this procedure we obtain a class of integrable models describing interacting 
quantum spins with non uniform long range interaction. 

The quasi-classical expansion of the transfer matrix reads 

t(u) = f (0) + 77 f (1) (u) + rf f (2) («) + O^ 3 ) (21) 

The property [£(«), t(t>)] = ensures the existence of hierarchy of integrable 
models in the quasi-classical expansion. We have indeed 

00 

[i(u),i(v)\ = j2v l Ci(u,v) = o, 

1=0 



Lj(u) 



r(u) = 



cosh(pw) 
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which implies C[(u,v) = 0. We give the first five terms 



C {u, 


v) 


=[fW(u) 


,f(°)(t)) 


5 












d(u, 


v) 


=[f<°>(u) 


,f«(t;) 


+ 


[f (1) (w) 


,T<P>(V) 








C 2 {u. 


v) 


=[f<°>(u) 




+ 


ifW(u) 




+ 


[f (1) ( M ) 




C 3 (u. 


v) 


=[f<°>(u) 




+ 


[f^(u) 


,T®(V) 


+ 


[f (1) ( M ) 




C 4 (u. 


v) 


=[f<°>(u) 




+ 


[f^(u) 


,t®(v) 


+ 










[f (1) («) 




+ 


[f^(u) 




+ 







i(2), 



(22) 



From the expansion above, one finds that the first non-trivial term 
(i.e. which is not just a (D-number or an invariant of the algebra) gives rise to a 
family of commuting operators. For example, if = C— number (as usually is 
the case) the first class of integrable models is generated by [f^(u), f^(v)] = 
0. In the presence of boundary conditions corresponding to generic choice of £± 
it turns out that f W (u) is non-trivial. This is not what we wish, since f W (u) 
only contains non-interacting spins.Q In the next section we will see how the 
parameters £± can be suitably chosen to yield "trivial" fW(u), such that the 
first non-trivial term in the quasi-classical expansion of t will be f^(u), which 
yields a spin-spin interaction. 

The r) expansion of the Lax operators is 



Lj(u) 



sinh(pu) 
P 



l + V lf\u)+ri 2 l) z, (u) 



2/(2), 



(23) 



where 



lf\u) 



lf\u) 



pcoth(pu)S?(T + , - — - 
J sinn ypu) 



(24) 



1 In general, the n-th order terms will contain up to n-body terms. 
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Thus, the monodromy matrix is: 



T{u) 



i 2 2 
+ P f] 



E 



[ ~ sinn [puj ) v J j j 

cosh (pu~) cosh (pu^)SjS^Jl + Sj~S^<r~a + + SjS^<T + cr^ 
sinh (pu~) sinh (pu^T ) 



+ 



det 9 T(-u 4 



^ E (^j) ^ + irrelevant terms j , 
r]/2)~P 2 (—u) 1 — jo?? coth (p»|) 



cosh (jm^)S z j <r x +Sf<r-+SjtT* 



i 2 2 
+ P 1] 



E 



cosh (put) cosh (put)S^S^ + S+S^ct-(t + + SjS^a+tr- 
sinh (ptij") sinh (pujj") 



+ 



1 ~ 2 

+ 2 E (^j) ^ + irrelevant terms 

3 



where we defined uf = u ± Zj, P(u) = FIj sinh(pn^) /p. The irrelevant terms 
are either (D-numbers or off-diagonal matrices, contributing to the trace with 
terms order r/ 3 . The expansion of K± reads 



1 ^sinh[p(M + d 0) )] 
~P\ -sinh[p( M -d 0) )], 

W+d 1) )«»h|p(«+d 0) )] 



+ 







(^-dVosh^-eni; 



(0)> 



where 5± + is the usual Kronecker 5, and we took into account the expansion 



(25) 



We define for convenience 2£ = ^ 2S = +^ 0) , and 2q = ^ 

Then the terms in the expansion of the transfer matrix given in Eq. (21) read 



11 



f (0) («) = ^P(u)P-\-u) (cosh (2pu) cosh (2pE) - cosh (2p-)) , 

f (1) (u) = C-numbers + ^P(u)P _1 (-u) x 
x sinh (2pw) sinh(2pE) ^(coth(pwj) + cot h (pit t)^? , 

3 

f (2 \u) = C-numbers + P(w)P _1 (-w)< 
2q sinh (2pu) cosh(2pE) ^ ^coth(pwj) + coth(pw+) J Sj + 
+ ^cosh (2pu) sinh (2pE) - sinh (2p£)^j ^ ^coth(pu~) + coth(pw+)^ 
+ - sinh (2pu) sinh(2pE) coth (pulf) coth (pu^)Sj + 

jr'fc 

+ - ^cosh (2pu) cosh (2pE) — cosh (2p£)^j x 



jk 



+ 



coih(pUj ) + coth(pw+)j ^coth(pw fc ) + coth(pn^) j SjS%+ 
2 (sinh(p-u~) sinh(p-u^) sinh(p-u^) sinh(p-u^)) ^ j ^ k ^ k ) 



2 

- if cosh (2pw) cosh (2p£)-cosh (2pE) ) ^ — - : : ' ; 



+ 



2 V 



~£ sinh (ptij- ) sinh (pn^ ) 



-j- ^ ' h (2 ) ' h (2 <£) k j 

2 * Z_/ gjjjjj (pit t) sinh (pu^) 



(26) 



As can be seen by Eqs. (22) , the operators (u) commute with each other if 



i(3)/ 



i(3)/ 



i(l), 



. 



A sufficient condition for this relation to be fulfilled is that f^(u) is just a 
(D-number. This requires that E = 0. 



3.1 Classical boundary 



At first, we assume classical boundary. With this term we mean that the 
boundary parameters £± are assumed to be independent of rj, i.e. ~ — 0. This 
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case was analyzed by one of the authors in Ref . [20] . The commuting family of 
operators f^(u) given above reduces to 



f(S 2 V) = 2P(u)P-\-u) sinh 2 (2pu)\J2 

1 





cosh (2pu) — cosh (2pzj) 3 



+ E 



i 

2 
1 



~k (cosh(2pw) — cosh(2pZj)) (cosh(2pu) — cosh.(2pz k )^ 
cosh(2p2j) + cosh(2pz fc ) - 2 cosh(2p£) J S*S k + 



+ - \cosh\p(zj + z k )\ - cosh(2p£) cosh[p(zj - z fe 
- 1 sinh(2p0 sinh[p(^ - z fc )] (fij"^ - 



)]) (5+5 fc - + SJS+) + 



where we dropped the (D-numbers and the Casimir coming from the term j — k 
in the sums. A finite subset of w-independent operators in involution can be 
obtained taking the limits u — > Zj of f( 2 \u), and dividing by the factor 



smh(2pz j )P 1 {-z j ) \\ sinh[p(^ - z k )\/p . 
Mi 



They are 



f 0j = J2 S k S j + 



■jf cosh(2pZj) — cosh(2pz fc ) 



cosh(2pZj) + cosh(2j9z fc ) - 2cosh(2p£) )S-S k + 



+ (cosh[p( Zj + z k )\ - cosh(2p0 cosh[p(^ - z fc )]) (5+4" + SjSt) + 
- sinh(2 P sinh[p(^ - z k )\ (S+ S k - SjSt 



(27) 



The fj form a complete set, in the sense that any f^(u) can be built from 
them according to the formula 



f (2) (w) = 2 P(u)P-\-u) sinh 2 (2pu) ^ 



1 



cosh(2pu) — cosh.{2pZj] 



TQj 



We notice the term J2kS k Sj = S" 2 ^ in the operators fj. It describes a self- 
interaction of the spins with the magnetic field generated by the spins them- 
selves. In the next section we shall see how to add an external magnetic field. 
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The eigenstates of operators (27) are given by 



M 



m=US-(e a )\H) , 

a=l 



(28) 



where 



S-(u) = £ + y + k d 

^ cosh(2pu) — cosh(2p;Zj) J a?7 



r;=0 



The rapidities e a fulfill the first order term in the expansion of Eqs. (19) 
around i] = 



y L - 

cosh(2pe Q ) - cosh(2pe^) 



E 



1/2 



(29) 



^ cosh(2pe a ) — cosh(2pzj) cosh(2pe a ) — cosh(2p£) 



. 



Putting cosh(2pe a ) = exp(2_E a ) + cosh(2p£) and cosh(2p2j) = exp(2u>j) + 
cosh(2p£), the equations above reduce to the modified Gaudin's equations 
presented in Refs. [15]: 



coth(£ a - E p ) - S j coth(E a - W j) + S z = 



(30) 



The eigenvalues are 



r 0j = Sj \ S z + y S k coth[(wj - w k )] - y c °th [(wj - E a ) 



(31) 



3.2 Non-classical boundary 



In this section we show how to include a scalable term proportional to Sj in the 
operators fj. Such a term is crucial for physical applications since, as we shall 
see, it allows to introduce a non-uniform magnetic field in the Hamiltonian. 
Furthermore, when the spins are realized by pairs of time-reversed electrons 
iS? = — + fiji — 1), Sj~ = Cj^cji a non-uniform magnetic field corresponds 
to a kinetic energy term (see section 4). In order to reach our goal, we have 
exploited the fact that £± can depend on 77, i.e. + ^ is not necessarily 
zero. We refer to this kind of boundary conditions as a non-classical boundary. 
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Thus, we put q ^ 0. We obtain 

f (2 \u) = 2P{u)p-\-u) sinh 2 (2pu){J2 

+ E 



2q + S z 



cosh (2pu) — cosh (2pzj) 3 



S*+ 



i,k (cosh(2p-u) — cosh(2pzj)) (cosh(2pu) — cosh(2pzk 
i ^cosh(2p^) + cosh(2pz fe ) - 2 cosh(2p£)^ S]S%+ 
+ i ^cosh[p(^ + z k )\ - cosh(2p0 cosh[p(^ - z fc )]^ + 5-^+) + 

- \ sinh(2 P sinh[p(^ - * fc )] - 5"^) | . (32) 

The integrals of motion are obtained again taking the limits u — > Zj, dividing 
now also by 2q + S z : 



1 



cosh(2pZj) — cosh(2p^ fc ) 



cosh(2pZj) + cosh(2pzfc) — 2 cosii^tj )>~>ju k 



h(2pS) y jS?SZ+ 



+ (cosh[p(^ + z k )\ - cosh(2p0 cosh[p(^- - ^)] J + + 



- sinh(2 P sinh[p(^- - - 



(33) 



where we put J(S Z ) = J/(l - JS 2 ), with J = -1/(2?). The operators t^ 2 \u) 
can be built from the fj according to 

f( 2 \u) = 2(2? + S^P^P'^-u) sinh 2 (2pu) V — - - . 1 — f ? . 

w v 7 v 7 v 7 i cosh(2pw) - cosh(2p^) J 

For real J, the r,- are Hermitian if Zj are real and £ is pure imaginary, or 
vice-versa. Their eigenstates are still given by Eq. (28) 



M 



m = n s '(^) \ H ) > 

Q = l 



(34) 



where 



{U> Y ~ cosh(2pw) - cosh(2p^) " j K dr? 



??=o 
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The difference with the previous subsection, is that the first order term in the 
expansion of Eqs. (19) around 77 = contains an additional term 



y L - 

cosh(2pe Q ) - cosh(2pe /3 ) 



-E 



+ 



i(i + 1/ j) 



(35) 



cosh(2pe a ) — cosh(2pzj) cosh(2pe a ) — cosh(2p£) 



Putting cosh(2pe a ) = exp(2E a ) + cosh(2p£) and cosh(2pzj) = exp(2vjj) + 
cosh(2p£), the equations above reduce as well to the modified Gaudin's equa- 
tions presented in Refs. [15]: 

coth(£ a -E p )-Y, Sj coth(£ Q - Wj ) + j^r = . (36) 



The eigenvalues are 



T j — Sj 



1 - J{S Z ) [ ]T Sk coth[( Wj - w k )} - ]T coth [( Wj - E a )\ 

k^j a 



(37) 



There are parameterizations yielding rational Bethe equations. Among these, 
we consider 

Y~ — - = cosh(2pzj) — cosh(2p£) , ^ ^ = cos h(2pe a ) — cosh(2j9^) . (38) 

J CI 

Then the eigenvalues are 



Tj = Sj[i- j(s z ) XjXk + J(s z ) 1 XjXa 



The Bethe equations (35) become 
1 



Xj Xk 



E 



E 83 



j A a Xj 



1 /1 + J(S*)(1 + S Z ) 1- J(S Z )(1 + S Z 







(39) 



(40) 



2 J(S Z ) \ l + \ a l-X a 
In this form, they admit a two dimensional electrostatic interpretation [26]. 



3.3 On the equivalence with Gaudin's model in external magnetic field 



We show that, at £ = and p^ = in/2, these integrals of motion are equivalent 
to the modified Gaudin's Hamiltonians introduced in Ref.[15]. The integrals 
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of motion reduce to 



a -qz _ KQ^sr) cosh ( 2 ^i) + cosh(2p* fc ) T 2 dz dz 
j j [ b) ££\ cosh(2^) - cosh(2pz fc ) ^" fe + 

cosh[p(z j + z fc )] =F cosh[p(^ - z fc )] 



cosh(2p2j) — cos 



h[p{ Zj - z k )\ ( q+q -, q - q +\\ 



where the upper sign refers to £ = and the lower one to p^ = in/ 2. We make 
the change of variable sinh (pzj) = expwj, if £ = 0, cosh(pzj) = expwj, if 
p£ = in/2, obtaining 

f, = s*-j(s z ) gjoothK- w k )s;s' k + 2 ginh( ^ _ wk) (s+s k ~ + s 7 st) 



Thus, apart a sector-dependent rescaling of the coupling J — > J(S Z ) — J (1 — 
JS* 2 ) -1 , the operators given above are equivalent to modified Gaudin's Hamil- 
tonians. 



3.4 Construction of the Hamiltonian 

The simplest Hamiltonian that is possible to be built up is a first degree 
polynomial in fj with arbitrary real parameters hj 

H= ^ 2H > f * = % ~ J(SZ) ^ cosh(2p^) - cosh(2p* fc ) 

(cosh(2pz j ) + cosh(2pz k ) - 2 cos [2pt]) S*S z k + 
+ (cosh\p(zj + z k )]-cos(2pt) cosh\p( Zj - z k )]j + 

- i sin(2pt) smh\p(zj - z k )\ (S+S k - S~S^ 
where we put £ = it, with real t. 



4 Second quantized Hamiltonians 



In this section we employ the fermionic realization of su(2) to write the inte- 
grable models we found Eq. (42) in second quantization. The two orthogonal 
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Dj and /^'-dimensional realizations are 

Da , D 



* - f 1 

^ + = E c ^mt> *7 = , = ^ E ( x - *W - *W) , 

(43) 



9 

<5, =1 ^ 5,=1 



and 



■ - - \ t 1 ' 



S^E^w , ST = (s+y , 5f = - E (n IpiT - n Ipii ) , (44) 



pi=i 2 ft=i 

where operators c, c^, and n = c^c are fermionic operators. We arbitrarily 
grouped the levels in the subsets j = 1, . . . , Qk, each containing Dj levels, and 
in the subsets 1 = 1,... ,fls, each containing D[ levels. The maximum values 
of the z components of the spin are Kj = Dj/2 and Si = D[/2 respectively. 
Thus, a level a will be characterized alternatively by the pairs (J (a), 5j( a )(a)) 
or (1(a), pi( a \{a)). We write a Hamiltonian of the form 

H = H K + H S + E Q , (45) 

where E is a constant and 

H K = -Y: 2 Vn (K) + £ 9jj k) , (46) 

3=1 3=1 

H s =Y J 2C l r l (S)+j:jrS 2 l , (47) 
i=i i=i 

where operators r(0) , O = K,S are defined in Eq. (33). Due to the orthogo- 
nality of the realizations (43), (44) we observe that [H K , H s ] = 0. Furthermore, 
Hk and Hs are block-diagonal, and their common eigenstates are the direct 
product of the eigenstates of Hk and of Hs, each restricted to the subspace 
corresponding to one of its blocks [16]. The integrability together with the 
exact solution of the Hamiltonian (45) follows from the integrability of each 
H K , H s proved in Section 3.2 and from Eqs. (34), (36), and (37). 

Finally, the second quantized form of the Hamiltonian (45) reads 

U ab (n a1 + n al )(n b ^ + n H ) + g ab c^cl l c H c b ^+ 



H = E £ aafl a a + E 
au ab 



(48) 

+ J ab(^ ~ n al )(n b] - n H ) + J^c^cJ^c^ 



where Q number the levels, a, b — 1, . . . , Q and c a>a = Cj^ ; s j(a) ( a ),a] the constant 
in Eq. (45) turns out to be E = J2j Dj Ej + J2jk DjD k Ujk- 
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The kinetic energy term reads 



aa 



We choose a partition — in equivalence classes — of the single particle levels 
in such a way that all levels having the same value of e a = \{e a ^ + £ a j.) belong 
to the same class (hence we write Ej instead of e a , where j individuates the 
class) Analogously, a second partition is defined in such a way that all 
the levels having the same value of ( a = \{s a -\ — £ jj belong to the same 
class (hence we write the common value as Ci) The couplings between 
levels a and b depend only on the equivalence classes of the two levels. For 
j = j(a) 7^ k = k{b) and I = 1(a) ^ m = m(b), they are 

9 r (iczsi n „ ^ cosh[p(zj + z k )}-cosh[p( Zl -z k - 2it K )} 
g ab = g jk = 2J K (K )( Vj - rj k ) cosh{2pZj) _ CO sh(2pz fe ) ' 

co$h(2pZj) + cosh(2pz fc ) — 2 cos (2ptx) 

cosh(2pZj) — cosh(2pz k ) 
cosh[p'(yj + y m )} -cos\i[p'(yi -y m + 2it s )} 
cosh(2p'yi) - cosh(2p'y m ) 
cosh(2p'yi) + cosh(2p'y m ) — 2 cos (2p'ts) 
cosh(2p'yi) - cosh(2p'y m ) 

(49) 



4U ab = 4U jk = J K (K z ){ Vj -r] k ) 

Jab = Jim = -2Js{S Z ){(l - Cm 
Jab = Jim = -Js{S Z )((l - (m) 



For j = A;, we have the relation g^ = 4:Ujj, and can be chosen arbitrarily.^ 
Analogously, for I = m, we have Jfi = Jf* 



5 Conclusions 



In this paper we have studied integrable disordered vertex models in presence 
of boundary reflections. The quasi-classical expansion of the models has been 
thoroughly investigated. This expansion produces a hierarchy of models which 
import the integrability of the original vertex models. The extraction of the 
energy from the generating functionals (which is unfeasible, in general) is 
very simplified by the quasi-classical limit and the Hamiltonian is given as 
polynomial of the integrals of motion. The class of models we obtain describes 
interacting spins with non uniform couplings, and in a non uniform external 
magnetic field. In this sense the present models generalize those ones found 

2 Vj = £ j + ^YLkDkUj k + gjj/2 have to be determined consistently; they must 
satisfy a system of linear equations, as discussed in [27]. 

3 The partitioning of the levels that we chose above guarantee that the interaction 
does not vanish between levels having the same value of e a or Ca- 

4 In particular, they can be chosen in such a way that r]j = Ej [28]. 
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in Ref. [20]. On the other hand, these Hamiltonians constitutes also a one- 
parameter (t in the text) extension of the class of models found in Ref. [15] 
that are recovered when the boundary terms give rise of twisted periodic 
boundary conditions. As a result, the integrability of these latter models has a 
firm ground within the Sklyanin procedure [19]. The presence of the external 
magnetic field is an effect of the boundary terms which are assumed, in turn, 
quasi-classical. We also obtained the exact solution of the class of models 
presented here through Algebraic Bethe Ansatz. The Bethe equations can be 
recast in a form which allows the electrostatic analogy as was done in Ref. [26]. 

An important point is that the models apply to any spin Sj (not only to spin 
1/2). The reason is that, in the present case, the integrals of motion can contain 
only spin S operators since the quantum algebra su q {2) reduces to su{2) in 
the quasi-classical limit (whatever the dimension of the representation is). 

By realizing the spin operators in terms of fermions, the class of models we 
found describes confined fermions in degenerate levels with pairing force in- 
teraction. 



A Inhomogeneous vertex models 

VM are models of 2D classical statistical mechanics. They consist in a (K x N) 
array of vertices (see Fig.A.l), where the nearest neighbours are linked by 
horizontal and vertical legs. The legs can be of several species, each identified 
by a number hij = 1, • • • ,-f/j, for the horizontal legs of the i-th row, and 
Vij = 1, • • • , Vj for the vertical ones of the j-th column (the number of species 
can depend on the row or column; in this case, we have an inhomogeneous 
model). Here we are using the convention that the pair individuates the 
vertical (horizontal) leg above (left of) the vertex A statistical weight 

1 2 3 N=4 



I 2 S> S> SI SI 

© ! © 

| K=3 (1 S> (I 

V l_J.ll 

Fig. A.l. (a) A vertex configuration, (b) Numbering of the lattice. 

w(legsi t j]i,j) = exp [— (3 e(legsij;i, j)] is assigned to each vertex, depending 
on the legs configurations (legsij = {hij, h i+ ij, v it j, i^+i}) around it. If the 
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weight depends explicitly on the position of the vertex, we have a disordered 
model. The goal is to find the partition function Z = J2{i egs } T\i,j w(legsij; i, j). 

To each vertex, one can associate a matrix, whose elements are the weights 
corresponding to the possible configurations of legs, in the following way: fix 
the horizontal legs around site to their minimal values, say hij = /li+ij = 
1; then vary the values of the vertical legs, associating the upper one to a 
row, and the lower one to a column; a (Vj x Vj) matrix is thus obtained, 
which we indicate by Lj' (i), whose entries are the weights corresponding to 
the possible legs configurations with horizontal legs fixed to 1; then repeat 
the procedure changing the values of the horizontal legs, associating the left 
one to a row, and the right one to a column; a block matrix Lj(i) is finally 
obtained, which is conventionally called the Lax operator. It is a (Hi x Hi) 
matrix whose entries L^' 1 ' 1+1,3 (i) are in turn (Vj x Vj) matrices, i.e. operators 
over the linear space Tij. The partition function of the (1 x N) lattice with 
periodic boundary conditions in vertical and horizontal direction is, in terms 
of the Lax operators, Z\ = TryTrf/{Li(l) • ■ ■ L N (1)} = Tryt(l), where by Tr H 
we mean the trace over the horizontal space, and by Try the trace over the 
vertical ones; we introduced the transfer matrix t(i) = Tr#{Li(i) • • ■ L^(i)}, 
where the hat is meant to remind that the transfer matrix is an operator over 
T~t = ®jT~Cj, the direct product of the linear spaces associated to the vertical 
legs. For a (K x N) lattice, the partition function is Z = Try{i(l) . . . t(K)}. 
If [t(i),t(i')} = Vi, i', it is possible to simultaneously diagonalize the t(i), 
obtaining Z = J2 r Ilf=i tr(i), where t r (i) is the r-th eigenvalue of i(i). In this 
case, the VM is exactly solvable. 

From the t(i), it is commonly possible to extract many-body Hamiltonians of 
interest Thus, a given exactly solvable vertex model corresponds uniquely 
to a family of commuting many-body operators. 

It turns out that the transfer matrices commute with each other, and thus 
the corresponding vertex model is exactly solvable, if Hi = Hi> = H, Vi, i', 
and a family of (H 2 x H 2 ) matrices, the i?-matrices, exists, such that the Lax 
operators obey the relation 

R(i,i') Lj(i) <g> Lj(if) = Lj(%) ® Lj(i) R(i,i') . (A.l) 

Given a i?-matrix, this is a very strict requirement, which in general implies 
that many legs configurations are not allowed, i.e. their weight is zero, while 
the allowed ones are related to each other by some parametrization. 

A relevant case is when the dimensions of vertical and horizontal space are 
equal, and they do not depend on the row or column: Hi = Vj = 2S + 1. 
Then, the i?-matrices are but the Lax operators where the matrix elements 
have been written down explicitly in their matrix representation. It turns out 

5 In general, such Hamiltonians are not by any means related to the Hamiltonian 
of the VM. 
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that the entries of the Lax operators are matrices belonging to the (2S + 1)- 
dimensional realization of su(2), i.e. spins over the Hilbert^] space Tij. There is 
the drawback that the ^-matrix is difficult to determine and to handle, since 
its dimension increases very fast with S [29,30]. A technique to build larger 
.R-matrices using the (4 x 4) i?-matrices (the simplest ones) as building blocks 
was devised by Kulish, Reshetikhin, and Sklyanin [31]. In the present paper, 
by means of the quasi-classical expansion, we will build up operators for spins 
higher than 1/2 still making use of (4 x 4) ^-matrices. 



B General K. matrix 



In this appendix we construct the Hamiltonian when the general solution of 
the reflection equation (15) is considered [33,34]. In this case we have 



K + (u) = - 
V 



sinh [p(tt + r] + £ + )] k + sinh[2p(tt + 77)] 
k + sinh[2p(w + 77)] — sinh[p(u + 77 — £+)] 



(B.l) 



Since we want f^ 1 ' to be a (D-number, we must impose k + ~ ir/c. Thus, the 
final effect of the general reflection results in an additional term in the second 
order of the transfer matrix 



f(2)( u ) ^ f (2 \u) + —P(u)P~\-u) sinh 2 (2pu)< 
P 

(B.2) 

y, smh.\p(zj - Q] g+ _ smh[p(zj + Q] ^_ 

j cosh(2pw) — cosh(2pZj) 3 ~" cosh(2p-u) — cosh(2pZj) ■* 



The Hamiltonian is again built according to 

H = J2 2h 3 r 3 (B.3) 
j 



6 It is a finite-dimensional vector space. We denote it as Hilbert space in foresight 
of its interpretation as a quantum space. 
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where the integrals of motion are 



7J = (1 - JS>)S; - icJ ^"h[Pfe-% + _ rinhlpfe+Ol g-^ + 



cosh(2p2j) — cosh(2p2;fc) 



(cosh(2p^) + cosh(2pz fc ) - 2cosh(2pf))S?S£+ 
+ (cosh[p(^ + * fc )] - cosh(2pO cosh[p(^ - z k )]) (S+S^ + + 

- sinh(2p£) smh\p(zj - z k )\ [SfS^ - Sj St) ■ 

(B.4) 

We point out that the Hamiltonian is hermitian for real c and £ = it with 
real t. The diagonalization of this class of Hamiltonians might be achieved by 
functional Bethe ansatz[6]. Nevertheless, it seems worth to study the models 
(B.3), (B.4) since their potential application to condensed matter (see also 
Eqs (43), (44)). 



C Generalization to su(n) 

In this appendix we briefly discuss on a generalization of the Gaudin model to 
the su{n) case. We can define the Gaudin model for other Lie algebras (see, 
e.g., Ref. [32] for recent works). Following the method depicted in section 3 we 
obtain a modified Gaudin Hamiltonian to include a scalable term proportional 
to the Cartan generators of su{n) (as far as we know, previously obtained 
models do not contain this term). 

The trigonometric i?-matrix for su(n) chains is given by 

in i n 

R(X) = -J2 sinh(p(A + r]))E aa ®E aa + -Y, sinh(pA)E aa ® E bb 

P a=l P a+b 

1 " 

+ -Y e - pXsgn{a - b) smh(pr])E ab <g> E ba , (C.l) 
P^b 

where E ab denotes n x n matrix with unity at (a, b) element. They satisfy 

[E ab , E cd ] = ^6 bc E ad — 5 da E cb ^j . 
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The corresponding diagonal solution of the reflection equation (14) is [33] 



K_(X) 
K + (X) 



- e pX sinh(p(£_ - \))E aa + - E e~ pA sinh(p(£_ + X))E aa , 

P a=l P a=i-+l 



~ E e- pX ~ 2pria sinh(p(£+ + \))E aa 

P a=l 

1 " 

+ - E e pX+Pv{n ~ 2a) smh(p(£ + - r]n - \))E a 

P a=i + +l 



(C.2) 



where t± is arbitrary, 1 < £± < n. Hereafter we set = for simplic- 

ity. With these if-matrices, the Hamiltonian of the su(n) homogenous spin 
chain with nearest neighbour interaction with open boundary was computed 
in Refs. [35,36] by using the formula (20). The Hamiltonian with the long 
range interaction is constructed following the procedure presented in section 
3.4, where the constants of motion are calculated by the formula (21). They 
read 



,=2E^r 



a=l 



+E 

a=l 



^ C oth(p^-z j ))+e+' coth(p^+ Zj ))+(n-£) 



sinh(2p Zj) 



sinh(p(£ - Zj )) sinh(p(£ + Zj )) 



+ E 

a=e+i I 



^ coth(p(£ + Zj )) + coth(p(e - Zj)) 



(1) 



n 

(coth(p( Z j -z k )) + coth(p(zj + z k ))) ]T E?E% 



+ E 

k+3 

n eP(zj - Zk)sga{a - b) I n sinh(p(e + %)) &a ^ b 

^ b smh{p{zj-z k )) 1 k ^ lb ^ +1 smh(p^-Zj))smh(p(zj + z k )) 3 k 

^ Z - 1 sinh(p(e + ^))sinh(p(2 ;i + ^)) 3 k 



a=i+l b 
e p(zj+z k )sga(b-a) 



n p(zj+z k )sgn(b-a) ^ 
_|_ \ " jjjbajjjab _|_ \ " _ jjjbaj^ab 

a j^ lS mh(p(zj + z k )) 3 k a ^ +1 smh(p(zj + z k )) 1 k 

a^b a^b 



(C.3) 



where E% b are site- A; su(n) operators. 



The spectrum of fj is given as limit of the eigenvalues of the sw(n)-transfer 
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matrix (Eq.(6) of [35]) and with u — > zf 



Tj = (coth\p( Zj + 0] - COth(pO) - (coth[p(^- + 0] + COth(pO) 



sinh(2p Zj) e 



-2zj 



+ (£- n) ) . J ' -rr + 1 - n 



sinh[p(zj — £)] sinh[p(,Zj + £)] sinh(2zj) 

jv , 

+ Y { coth[p(zj + z k )\ + coth[p(z,- - z fc )] 

Mi , 

+ £ coth[p(^ + ei 1} )] + coth^ - e^)] 



k 



(C.4) 



The Bethe ansatz equations can be obtained in the same limit of a result in 
Ref. [35] , and we have (for a — 1, 2, . . . , n — 1) 



2 £ (coth[p(ef + e«)] + coth[p(ef } - <£>)]) 
Mi v 7 

+ <S M (n + £<_ 1} - d 1} ) (coth[p(z - 0] + coth[p(z + 0]) 

Ma+l , v 

= £ (coth[ P (ef } + et +1) )} + coth[p(ef - e^)]) 

M a _i 

+ 2 (coth[p(ef + et 1] )] + coth[p(ef } - e^)]) (C.5) 
Here we assume ef^ = Zj and M = N, M n = 0. 

The fermionic models can be obtained by using the fermionic realization 

Ef = c\ a c hh -U ah , (C.6) 

with a constraint 

^ C jr':« c i' a = ^ ' 



a 
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